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ABSTRACT 

Astrometric  measurements  are  presented  for  seven  nearby  stars  with  previously  detected  planets:  six  M  dwarfs  (GJ 
317,  GJ  667C,  GJ  581,  GJ  849,  GJ  876,  and  GJ  1214)  and  one  K  dwarf  (BD-10  -3166).  Measurements  are  also 
presented  for  six  additional  nearby  M  dwarfs  without  known  planets,  but  which  are  more  favorable  to  astrometric 
detections  of  low  mass  companions,  as  well  as  three  binary  systems  for  which  we  provide  astrometric  orbit 
solutions.  Observations  have  baselines  of  3  to  13  years,  and  were  made  as  part  of  the  RECONS  long-term 
astrometry  and  photometry  program  at  the  CTIO/SMARTS  0.9  m  telescope.  We  provide  trigonometric  parallaxes 
and  proper  motions  for  all  16  systems,  and  perform  an  extensive  analysis  of  the  astrometric  residuals  to  determine 
the  minimum  detectable  companion  mass  for  the  12  M  dwarfs  not  having  close  stellar  secondaries.  For  the  six  M 
dwarfs  with  known  planets,  we  are  not  sensitive  to  planets,  but  can  rule  out  the  presence  of  all  but  the  least  massive 
brown  dwarfs  at  periods  of  2-12  years.  For  the  six  more  astrometrically  favorable  M  dwarfs,  we  conclude  that  none 
have  brown  dwarf  companions,  and  are  sensitive  to  companions  with  masses  as  low  as  1  Mjup  for  periods  longer 
than  two  years.  In  particular,  we  conclude  that  Proxima  Centauri  has  no  Jovian  companions  at  orbital  periods  of 
2-12  years.  These  results  complement  previously  published  M  dwarf  planet  occurrence  rates  by  providing 
astrometrically  determined  upper  mass  limits  on  potential  super-Jupiter  companions  at  orbits  of  two  years  and 
longer.  As  part  of  a  continuing  survey,  these  results  are  consistent  with  the  paucity  of  super-Jupiter  and  brown 
dwarf  companions  we  find  among  the  over  250  red  dwarfs  within  25  pc  observed  longer  than  five  years  in  our 
astrometric  program. 
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1.  INTRODUCTION 

During  the  brief  history  of  extrasolar  planet  investigations, 
our  understanding  of  the  relative  populations  of  different  types 
of  planets  has  been  limited  by  the  observational  biases  of  the 
techniques  employed.  With  the  advent  of  sophisticated  transit 
searches  and  hypersensitive  radial  velocity  measurements, 
significant  progress  has  been  made  discovering  various  types 
of  planets  that  orbit  stars  with  periods  up  to  a  few  years.  Less 
progress  has  been  made  in  discovering  planets  in  longer  orbits, 
and  particularly  around  nearby  M  dwarfs,  which  account  for  at 
least  74%  of  the  stellar  population  within  10  pc  (Henry 
et  al.  2006).  M  dwarfs  offer  fertile  ground  for  companion 
searches,  as  Dressing  &  Charbonneau  (2013)  have  inferred  that 
a  high  fraction  of  M  dwarfs  host  terrestrial  planets  at  short 
orbital  periods.  Less  is  known  about  the  populations  of  Jupiter- 
mass  planets  and  brown  dwarfs  around  M  dwarfs,  particularly 
at  orbital  periods  longer  than  a  few  years. 

To  understand  how  M  dwarf  planetary  systems  form  and 
evolve,  we  must  probe  the  full  regime  of  companion  masses 
and  orbital  periods.  Transit  techniques  are  geometrically  biased 
toward  companions  with  small  orbits,  while  radial  velocity 
techniques  are  biased  toward  massive  companions  with  short 
periods  that  exert  large  gravitational  accelerations  on  their  host 
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stars.  Direct  imaging  techniques  are  limited  to  young,  giant 
planets  at  large  separations.  Astrometric  techniques,  which 
measure  the  positions  of  stars  on  the  plane  of  the  sky,  are  most 
sensitive  to  Jovian-type  planets  in  Jovian-type  orbits.  While 
radial  velocity  observing  programs  are  now  becoming  sensitive 
to  such  companions  (Bonfils  et  al.  2013a;  Montet  et  al.  2014), 
the  astrometric  results  presented  here  have  longer  observational 
baselines,  of  up  to  13  years.  Furthermore,  astrometry  can  detect 
companions  with  a  large  range  of  inclinations  and  orientations, 
and  allow  for  the  determination  of  orbit  inclinations  and 
accurate  companion  masses. 

To  date  the  majority  of  nearby  extrasolar  planets  around  M 
dwarfs  have  been  discovered  by  radial  velocity  searches,  which 
tend  to  select  the  brightest  M  dwarfs.  As  discussed  in  more 
detail  in  Section  3,  in  ground-based  imaging  programs  the 
brightest  targets  generally  have  the  noisiest  astrometric 
residuals  due  to  the  short  exposures  required  and  the  lack  of 
comparably  bright  reference  stars.  With  the  exception  of  GJ 
1214,  five  M  dwarfs  in  our  sample  were  found  to  have  planets 
using  radial  velocity  techniques,  and  are  among  the  brightest 
targets  in  our  astrometric  program.  An  extreme  case  is  the  K 
dwarf  BD-10  3166,  for  which  we  are  not  sensitive  to  sub-stellar 
companions,  but  for  which  we  provide  the  first  accurate 
parallax.  For  comparison,  we  have  included  six  additional  M 
dwarfs  that  are  less  bright,  less  massive,  and  closer,  and 
therefore  more  favorable  to  companion  detection  via  astro¬ 
metry.  To  calibrate  our  analysis,  we  have  also  included  three 
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confirmed  stellar  binaries  with  clear  photocentric  perturbations 
for  which  we  have  characterized  the  orbits.  These  binaries  were 
chosen  from  the  roughly  two  dozen  binaries  in  our  observing 
program  with  clear  astrometric  perturbations  because  we  have 
observed  multiple  orbital  periods,  and  can  most  accurately 
characterize  the  orbits. 

Astrometric  solutions  for  proper  motion  and  parallax  are 
given  for  each  of  the  16  systems  targeted,  plus  orbital  solutions 
for  three  binaries.  A  detailed  analysis  of  the  astrometric 
residuals  is  given  to  search  for  companions  to  the  12  M  dwarf 
systems  without  close  stellar  companions.  Periodograms  of  the 
astrometric  residuals  have  been  generated,  along  with  detection 
limits  based  on  simulations  of  10  million  hypothetical 
companions  to  each  star.  These  are  the  first  results  of  a  larger 
RECONS  survey  for  companions  orbiting  more  than  250  red 
dwarfs  within  25  pc  for  which  we  have  at  least  five  years  of 
coverage.  As  observations  continue,  this  sample  will  grow, 
further  constraining  the  population  of  brown  dwarf  and  super- 
Jupiter  companions  in  long  period  orbits  around  M  dwarfs. 
Finally,  to  provide  context  for  these  results  we  provide  a 
comprehensive  list  of  the  1 7  M  dwarfs  within  25  pc  having 
exoplanets  as  of  2014  July  1,  including  the  six  targeted  in 
this  work. 

2.  OBSERVATIONS  AND  REDUCTIONS 
2.1.  Astrometry 

The  0.9  m  telescope  at  CTIO  is  equipped  with  a  2048  x  2048 
Tektronix  CCD  camera  with  0^401  pixel-1  plate  scale  (Jao 
et  al.  2003).  Only  the  center  quarter  of  the  chip  is  used  for 
astrometric  and  photometric  observations,  yielding  a  6)8  square 
field  of  view.  Astrometric  exposures  are  taken  through  one  of 
four  filters,  Vj  (old),  Vj  (new),  RKc.  or  Jkc* 9  (hereafter  without 
subscripts,  and  the  V  filters  combined).  Depending  on  the 
brightnesses  of  the  science  targets,  reference  stars,  and  sky 
conditions,  exposure  times  vary  from  20  to  1200  s  for  targets 
with  9  <  VRI  <  19.  For  optimal  centroiding,  exposure  times 
are  set  so  that  either  science  or  reference  stars  have  maximum 
peak  ADU  of  ~50,000  (digital  saturation  occurs  at  65,537 
ADU).  Observations  are  almost  always  made  within  ±30 
minutes  of  a  science  target’s  transit  to  minimize  the  corrections 
required  for  differential  color  refraction,  as  described  in  Jao 
et  al.  (2005).  Three  to  five  frames  are  typically  taken  each 
night,  depending  primarily  on  the  exposure  time  required.  To 
enable  routine  calibration  of  the  science  images,  bias  and  dome 
flat  frames  are  taken  nightly. 

Instrument  setups  for  most  stars  have  been  kept  constant 
during  the  13  years  of  observations.  However,  we  have  used 
two  V  filters,  dubbed  the  “old”  Tek#2  V  filter  (Acentrai 
=  5438  A,  AA  =  1026  A)  and  “new”  Tek#l  V  filter  (Acentral 
=  5475  A,  AA  =  1000  A),  because  the  “old”  filter  cracked  in 
2005  February.  The  “new”  V  filter  was  used  between  2005  and 
2009.  The  “old”  V  filter  was  reinstated  in  2009  July  after 
confirming  that  the  crack  in  the  corner  did  not  significantly 
affect  astrometric  residuals.  As  discussed  in  Subasavage  et  al. 
(2009),  a  reliable  parallax  can  be  obtained  using  data  from  both 
filters  as  long  as  at  least  1-2  years  of  data  (depending  on 
observing  frequency)  have  been  taken  in  each  filter.  Reductions 
containing  both  “old”  and  “new”  V  frames  can  exhibit  offsets 
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of  a  several  milliarcseconds  (mas)  in  residuals  on  both  axes. 
This  has  been  mitigated  by  choosing  close-in  reference  stars, 
and  only  using  frames  taken  near  the  meridian.  In  total,  7  of  the 
16  systems  discussed  in  this  paper  were  observed  astrome- 
trically  in  the  V  filter.  Further  details  about  the  filters  and  their 
effects  on  the  astrometry  can  be  found  in  Subasavage  et  al. 
(2009)  and  Riedel  et  al.  (2010). 

The  paths  traced  on  the  sky  by  science  stars  result  from  the 
combinations  of  proper  motions  and  parallactic  shifts.  Details 
of  the  data  reduction  process  used  to  separate  these  motions  are 
given  in  Jao  et  al.  (2005)  and  Henry  et  al.  (2006).  Briefly,  we 
(1)  use  SExtractor  (Bertin  &  Arnouts  1996)  to  measure 
centroids,  (2)  define  a  six-constant  plate  model  to  find  plate 
constants  (given  in  Equation  (4)  of  Jao  et  al.  2005),  (3)  assume 
that  ensembles  of  reference  stars  have  zero  mean  parallax  and 
proper  motion,  (4)  solve  least-square  equations  for  multi-epoch 
images  (given  in  Equation  (5)  of  Jao  et  al.  2005),  and  (5) 
convert  from  relative  parallax  to  absolute  parallax  by  estimat¬ 
ing  the  distances  of  the  reference  stars  photometrically.  Our 
typical  centering  precision  is  2. 1-3. 5  mas,  or  0.5%-0.9%  of  a 
pixel,  depending  on  the  filter,  with  I  being  the  best  and  V  being 
the  worst.  To  correct  the  relative  parallax  to  an  absolute 
parallax,  photometric  distances  are  estimated  by  comparing 
VRI  colors  to  Mv  for  single,  main-sequence  stars  in  the 
RECONS  10 pc  sample  (Henry  et  al.  1997,  2006).  A  distance 
is  estimated  for  each  reference  star,  and  the  correction  to 
absolute  parallax  is  then  computed  using  the  weighted  mean 
distance  of  the  entire  reference  field.  The  uncertainty  on  the 
correction  is  determined  using  Equation  (6)  in  Jao 
et  al.  (2005). 

In  the  case  of  a  binary  with  a  given  combination  of 
magnitude  and  mass  differences,10  we  detect  its  photocenter 
orbit  around  its  barycenter,  in  addition  to  the  motions  due  to 
parallax  and  proper  motion.  Hence,  the  residuals  of  our  typical 
binary  star’s  parallax  reduction  are  significantly  offset  from 
zero.  In  order  to  get  a  better  parallax  result  and  calculate  the 
photocenter’s  orbital  elements,  we  first  carry  out  a  standard 
reduction  for  proper  motion  and  parallax.  We  then  fit  a 
photocentric  orbit  to  the  residuals,  i.e,  we  treat  these  residuals 
as  a  binary  orbit,  using  the  techniques  described  in  Hartkopf 
et  al.  (1989).  Based  on  the  orbital  elements  we  calculate,  this 
photocentric  orbit  in  components  of  right  ascension  (R.A.)  and 
declination  (deck)  is  then  removed  from  the  centroids  of  the 
science  star  at  each  epoch.  Finally,  we  re-calculate  proper 
motion  and  parallax  using  these  corrected  centroids.  The  final 
“cleaned”  proper  motions  and  parallaxes  are  the  values  given 
for  the  binaries  in  Table  1.  After  one  iteration,  the  parallax 
errors  are  reduced  to  those  typical  of  similar  program  stars,  and 
the  residuals  are  significantly  reduced  and  consistent  with  the 
mean  errors  found  for  our  overall  program. 

2.2.  Photometry 

VRI  photometry  was  obtained  at  the  CTIO  0.9  m  using  the 
same  instrumental  setup  used  for  the  astrometry  frames.  As  for 
astrometry  observations,  bias  and  dome  flat  frames  are  taken 
nightly  for  basic  image  calibration.  All  science  stars  were 
observed  at  airmass  <1.8.  Exposure  times  were  chosen  to  reach 


As  described  in  van  de  Kamp  (1967),  o  =  (71  -  3)a,  where  a  is  the 

fihotocentric  semimajor  axis  of  the  orbit  of  the  primary,  B  is  the  fractional  mass 
Mb/(Ma  +  Afg)),  0  is  the  relative  luminosity  (1/(1  +  10(0'4,Am)),  and  a  is 
the  semimajor  axis  of  the  relative  orbit  of  the  two  components.  The 
perturbation  we  detect  here  is  a. 
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Table  1 

Astrometric  Results 


Name 

a) 

R.A. 

(J2000.0) 

(2) 

Decl. 

(J2000.0) 

(3) 

Fil. 

(4) 

N  sea 

(5) 

N  frm 
(6) 

Coverage 

(V) 

Years 

(8) 

N  ref 

(9) 

n  (rel) 
(mas) 
(10) 

n 

(corr) 

(mas) 

(ii) 

n  (abs) 
(mas) 

(12) 

/i 

(mas  yr~') 
(13) 

P.A. 

(deg) 

(14) 

Vta„ 

(km  s-1) 
(15) 

Obs. 

Prec. 

(mas) 

(16) 

Perc. 

Cmpn. 

(17) 

Notes 

(18) 

Extrasolar  Planet  Hosts 

GJ  317 

08  40  59.21 

-23  27  22.6 

R 

5c 

75 

2009.04-2013.38 

4.35 

7 

64.04 

1.50 

65.54 

930.7 

330.5 

65.5 

4.99 

±  1.45 

±  0.50 

±  1.53 

zb  1.1 

zb  0.13 

BD-10  3166 

10  58  28.79 

-10  46  13.4 

I 

7s 

71 

2004.43-2011.50 

7.07 

6 

13.84 

1.50 

15.34 

185.9 

269.1 

52.4 

10.85 

±  3.04 

±  0.50 

±  3.08 

=b  1.5 

zb  0.67 

GJ  581 

15  19  26.83 

-07  43  20.1 

V 

14s 

267 

2000.58-2013.38 

12.80 

11 

157.67 

1.12 

158.79 

1224.3 

266.0 

36.5 

7.93 

96% 

a 

±  1.57 

±  0.17 

±  1.58 

:b  0.4 

zb  0.03 

GJ  1214 

17  15  18.92 

+04  57  50.1 

I 

4c 

80 

2010.39-2013.38 

3.00 

9 

68.20 

1.88 

70.08 

945.5 

142.0 

63.9 

5.02 

±  1.26 

zb  0.18 

zb  1.27 

zb  1.4 

zb  0.17 

GJ  667C 

17  18  58.82 

-34  59  48.6 

V 

1  Is 

140 

2003.52-2013.38 

9.86 

5 

139.38 

1.50 

140.88 

1154.1 

101.0 

38.8 

7.66 

93% 

zb  1.98 

±  0.50 

zb  2.04 

±  0.6 

zb  0.05 

GJ  849 

22  09  40.34 

-04  38  26.8 

V 

11s 

135 

2003.52-2013.39 

9.86 

5 

113.78 

2.27 

116.05 

1118.0 

90.8 

45.7 

9.53 

81% 

±  1.97 

zb  0.30 

±  1.99 

zb  0.5 

±  0.04 

GJ  876 

22  53  16.75 

-14  15  49.2 

V 

11s 

85 

2003.52-2013.39 

9.87 

6 

210.97 

2.14 

213.11 

1149.4 

125.7 

25.6 

8.43 

99% 

±  3.99 

±  0.57 

±  4.03 

zb  1.1 

zb  0.11 

Best 

Case  Targets 

GJ  1061 

03  35  59.72 

-44  30  45.5 

R 

13s 

194 

1999.62-2012.95 

13.32 

7 

269.92 

0.94 

270.86 

827.7 

117.7 

14.5 

7.59 

79% 

a 

±  1.29 

±  0.08 

±  1.29 

±  0.3 

zb  0.04 

LP  944-020 

03  39  35.25 

-35  25  43.8 

I 

8s 

59 

2003.95-2012.94 

8.99 

10 

154.53 

1.36 

155.89 

408.3 

48.5 

12.4 

2.13 

94% 

a 

±  1.03 

±  0.10 

±  1.03 

zb  0.3 

zb  0.07 

GJ  1128 

09  42  46.36 

-68  53  06.1 

V 

13s 

167 

2000.23-2013.12 

12.89 

8 

153.54 

0.73 

154.27 

1123.0 

356.1 

34.5 

3.17 

80% 

a 

±  0.75 

±  0.11 

±  0.76 

zb  0.2 

zb  0.02 

DENIS  J 1048-3956 

10  48  14.56 

-39  56  07.0 

I 

13s 

200 

2001.15-2013.27 

12.13 

11 

247.23 

0.85 

248.08 

1531.6 

229.5 

29.3 

2.92 

97% 

a 

±  0.60 

zb  0.10 

±  0.61 

zb  0.2 

zb  0.01 

SCR  1138-7721 

11  38  16.76 

-77  21  48.5 

I 

11s 

134 

2003.25-2013.27 

10.03 

12 

119.60 

0.81 

120.41 

2143.3  zb 

287.8  zb 

84.4 

4.20 

69% 

a 

zb  1.01 

zb  0.07 

zb  1.01 

0.4 

0.02 

Proxima  Cen 

14  29  43.02 

-62  40  46.7 

V 

14s 

205 

2000.57-2013.25 

12.68 

5 

766.41  ± 

1.72  ± 

768.13  ± 

3850.8  zb 

282.4  zb 

23.8 

4.83 

99% 

a 

0.91 

0.50 

1.04 

0.6 

0.02 

Confirmed  Binaries 

LHS  1582AB 

03  43  22.08 

-09  33  50.9 

R 

11s 

102 

2000.87-2012.94 

12.06 

7 

48.84  zb 

2.00  ± 

50.84  zb 

509.4  zb 

52.7  zb 

47.5 

3.88 

a 

1.18 

0.26 

1.21 

0.3 

0.06 

GJ  748 AB 

19  12  14.60 

+02  53  11.0 

V 

10s 

154 

2004.45-2013.39 

8.95 

11 

97.77  zb 

2.22  ± 

99.99  zb 

1857.8  dz 

107.4  zb 

88.1 

5.80 

1.15 

0.41 

1.22 

0.5 

0.02 

LHS  3738AB 

21  58  49.13 

-32  26  25.5 

R 

12s 

151 

1999.64-2012.81 

13.17 

12 

50.82  ± 

1.40  ± 

52.22  zb 

535.2  =b 

229.1  zb 

48.6 

2.50 

a 

1.01 

0.21 

1.03 

0.3 

0.06 

Note.  Nsea  indicates  the  number  of  seasons  observed,  where  3-6  months  of  observations  count  as  one  season,  for  seasons  having  more  than  three  images  taken.  The  letter  “c”  indicates  a  continuous  set  of  observations 
during  which  multiple  nights  of  data  were  taken  in  each  season,  whereas  an  “s”  indicates  scattered  observations  when  one  or  more  seasons  have  only  a  single  night  of  observations.  Generally,  “c”  observations  are  better, 
(a)  Target  has  one  or  more  parallaxes  previously  published  by  RECONS.  The  values  here  supersede  those  earlier  values.  Perc.  Cmpn.  indicates  the  percentage  of  eliminated  brown  dwarf  companions  to  the  extrasolar 
planet  hosts,  and  the  percentage  of  eliminated  planetary-mass  companions  to  the  best  case  targets. 
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a  signal-to-noise  ratio  (S/N)  >100  for  science  stars  in  each  of 
the  VRI  filters.  Combinations  of  fields  that  provided  10  or  more 
standard  stars  from  Landolt  (1992,  2007)  and/or  E-regions 
from  Graham  (1982)  were  observed  several  times  each  night  to 
derive  transformation  equations  and  extinction  curves.  Further 
details  of  photometric  data  reductions,  the  definition  of 
transformation  equations,  errors,  etc.,  can  be  found  in  Jao 
et  al.  (2005)  and  Winters  et  al.  (2011). 

2.3.  Spectroscopy 

Spectroscopic  observations  used  to  the  provide  spectral  types 
in  Table  2  were  made  at  the  CTIO  1.5  m  using  the  R-C 
spectrograph  and  Loral  1200  x  800  CCD  camera  between  2003 
and  2006.  Grating  #32  was  used  in  first  order  with  a  tilt  of  15°  1, 
and  observations  were  made  using  a  2"  slit.  The  order-blocking 
filter  OG570  was  utilized  to  provide  spectra  covering  the  range 
of  6000-9000  A  with  a  resolution  of  8.6  A.  For  calibration, 
bias  frames,  dome  flats,  and  sky  flats  were  taken  at  the 
beginning  of  each  night.  Further  details  regarding  the  1.5  m 
spectroscopy  program  and  associated  data  reduction,  including 
assignment  of  spectral  types,  can  be  found  in  Henry 
et  al.  (2004). 

3.  ANALYSIS 

Astrometric  residuals  represent  the  deviation  in  a  target’s 
measured  position  from  the  solution  for  proper  motion  and 
parallax,  given  for  all  16  target  stars  in  Table  1.  The  residuals 
for  each  star  are  plotted  in  Figure  1 .  Each  filled  circle  represents 
the  mean  of  typically  three  to  five  frames  taken  in  a  single 
night,  with  a  corresponding  estimate  of  the  nightly  mean  error. 
Open  circles  represent  nights  with  only  one  frame,  which  are 
included  in  the  parallax  measurement  but  are  not  included  in 
the  following  analysis  of  the  residuals.  The  three  binaries  (GJ 
748AB,  LHS  1582AB,  and  LHS  3738AB)  exhibit  large, 
periodic  perturbations  in  both  R.A.  and  deck,  indicating  the 
presence  of  a  companion  in  each  case.  Because  the  systems  are 
unresolved  in  our  images,  we  calculate  a  photocentric 
semimajor  axis,  as  previously  discussed.  The  remaining  targets 
generally  have  flat  residuals,  although  some  are  more  scattered 
(GJ  849)  than  others  (GJ  1128).  Targets  with  the  flattest 
residuals,  such  as  DENIS  J1048-3956,  generally  have  at  least  8 
reference  stars  that  closely  and  evenly  surround  the  target,  and 
are  bright  enough  to  have  a  least  1,000  peak  counts  when  the 
target  has  50,000.  Fainter  targets  tend  to  have  the  flattest 
residuals  because  they  have  a  larger  number  of  suitable 
reference  stars  and  require  longer  exposures,  which  smooth  out 
short-term  seeing  effects  that  produce  PSFs  with  poorly- 
defined  centroids. 

The  residuals  were  analyzed  for  the  presence  of  companions 
using  Lomb-Scargle  periodograms,  shown  in  Figure  2,  gener¬ 
ated  by  the  method  given  in  Zechmeister  &  Kiirster  (2009). 
This  method  weights  each  data  point  by  its  estimated  nightly 
mean  error,  and  produces  a  normalized  periodogram  for  the 
data  sequence  with  power  ranging  from  0  and  1.  We  only 
attempt  to  detect  companions  with  periods  greater  than  two 
years  but  less  than  the  length  of  the  observations,  and  have  only 
generated  periodograms  for  those  ranges.  The  data  were 
oversampled  in  increments  of  one  day  to  create  smooth 
periodograms.  As  Frescura  et  al.  (2008)  have  shown,  over- 
sampling  the  periodogram  does  not  significantly  increase  the 
risk  of  false  periodogram  peaks. 


As  the  data  are  irregularly  spaced,  it  was  not  possible  to  use 
an  analytical  formula  for  the  false  alarm  probabilities  (FAPs). 
Instead,  empirical  FAP  distribution  functions  were  calculated 
as  outlined  in  Frescura  et  al.  (2008)  by  generating  period¬ 
ograms  for  10,000  sets  of  random  noise.  Each  of  the  10,000 
data  sets  has  the  same  number  of  observations  as  the  actual 
data,  with  the  same  observation  times  and  errors  as  the  actual 
data.  The  randomized  values  of  each  point  within  a  given  data 
set  are  assigned  from  a  normal  distribution  with  a  standard 
deviation  representing  the  distribution  of  offsets  from  zero  of 
the  actual  data  points. 

Each  of  these  10,000  periodograms  was  then  sampled  in  a 
grid  of  periods  with  five  day  increments  between  two  years  and 
the  length  of  observations  to  find  the  highest  power  occurring 
in  each  periodogram.  The  FAP  for  an  observed  periodogram 
peak  of  a  given  power  is  equal  to  the  fraction  of  the  10,000  data 
sets  with  peaks  greater  than  or  equal  to  that  power.  The 
normalized  powers  corresponding  to  FAP  values  of  0.1,  0.01, 
and  0.001  are  shown  as  dashed  lines  in  Figure  2.  This 
corresponds  to  10%,  1%,  and  0.1%  probabilities  that  the  peak 
is  due  solely  to  random  noise  and  the  cadence  of  observations. 
As  a  calibration,  periodograms  and  FAP  functions  were  also 
generated  for  the  three  confirmed  stellar  binaries  for  which 
astrometric  perturbations  are  clearly  evident.  These  binaries 
have  significant  periodogram  peaks  (FAP  <  0.001)  in  both  the 
R.A.  and  deck  axes  at  periods  corresponding  to  their  orbital 
periods,  while  the  remaining  13  targets  have  no  significant 
peaks.  This  indicates  that  our  periodograms  are  sensitive  to 
astrometric  perturbations,  and  that  we  have  not  detected  any 
significant  periodicity  in  the  residuals  of  the  remaining  12  M 
dwarfs  and  one  K  dwarf  without  companions. 

Given  that  there  are  no  companions  evident  in  the  residuals 
of  the  non-binary  targets,  we  aim  to  establish  lower  limits  for 
the  companion  masses  and  periods  to  which  we  are  sensitive, 
and  would  have  detected  were  they  present.  For  each  target  we 
ran  a  simulation  of  10  million  hypothetical  companions 
orbiting  each  star  with  masses  chosen  randomly  from  a 
uniform  distribution  between  0.5  and  80  MJup,  and  periods 
chosen  randomly  from  a  uniform  distribution  between  2  years 
up  to  the  length  of  observations.  Although  we  have  observed 
perturbations  as  short  as  1 .2  years,  our  practical  lower  limit  is  2 
years.  All  geometric  parameters  (inclination,  eccentricity,  time 
of  periastron,  longitude  of  ascending  node,  and  longitude  of 
periastron)  were  assigned  randomly  from  a  uniform  distribu¬ 
tion  of  all  possible  values,  including  eccentricity,  which  was 
allowed  to  be  as  high  as  0.99. 

The  program  then  calculates  the  astrometric  perturbation  that 
each  simulated  companion  would  induce  on  the  primary,  using 
a  primary  mass  calculated  from  the  mass-luminosity  relations 
in  Henry  &  McCarthy  (1993)  (Equation  (5a,  b))  and  the 
revised  relation  for  the  lowest  mass  stars  in  Henry  et  al.  (1999). 
For  DEN  J1048-3956  and  LP  944-020,  which  are  too  faint  for 
the  relations,  we  assume  masses  of  0.08  M0.  We  assume  that 
the  companion  does  not  contribute  significantly  to  the  overall 
flux  in  each  system,  so  that  the  photocenter  of  the  system  is 
concentric  with  the  primary  star.  This  assumption  proves 
problematic  when  considering  brown  dwarf  companions  to  the 
two  latest  M  dwarfs,  DEN  J 1048-3956  and  LP  944-020. 
However,  the  trigonometric  parallaxes  of  these  targets  agree 
well  with  their  photometric  distance  estimates  (Columns  12 
and  13  of  Table  2),  and  we  conclude  that  neither  have  nearly 
equal  luminosity  companions. 
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Table  2 

Photometric  and  Spectroscopic  Results 


Name 

V 

R 

1 

Nights 

J 

H 

Ks 

Spectral 

Ref. 

Mass 

Trig.  dist. 

Phot.  dist. 

No.  of 

Notes 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

type 

(M0) 

(Pc) 

(Pc) 

relations 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

Extrasolar  Planet  Hosts 


GJ  317 

BD-10  3166 

12.01 

10.03 

10.84 

9.58 

9.37 

9.19 

3 

3 

7.934 

8.611 

7.321 

8.300 

7.028 

8.124 

M3. 5  V 

K3.0  V 

1 

2 

0.35 

0.85 

15.26 

65.19 

± 

± 

0.36 

13.64 

9.70 

± 

1.53 

12 

a 

GJ  581 

10.56 

9.44 

8.03 

3 

6.706 

6.095 

5.837 

M3.0  V 

2 

0.30 

6.30 

± 

0.06 

6.60 

± 

1.03 

12 

GJ  1214 

14.71 

13.27 

11.50 

3 

9.750 

9.094 

8.782 

M4.5  V 

1 

0.14 

14.27 

± 

0.24 

12.42 

± 

2.00 

12 

GJ  667 C 

10.34 

9.29 

8.09 

3 

6.848 

6.322 

6.036 

Ml. 5  V 

2 

0.36 

7.10 

± 

0.10 

9.41 

± 

1.49 

12 

GJ  849 

10.38 

9.27 

7.87 

3 

6.510 

5.899 

5.594 

M3.0  V 

2 

0.42 

8.62 

± 

0.15 

5.73 

± 

0.92 

12 

GJ  876 

10.18 

8.97 

7.40 

3 

5.934 

5.349  5.010 

Best  Case  Targets 

M3. 5  V 

2 

0.27 

4.69 

± 

0.09 

3.46 

± 

0.54 

12 

GJ  1061 

13.09 

11.45 

9.47 

6 

7.523 

7.015 

6.610 

M5.0  V 

3 

0.11 

3.69 

± 

0.02 

3.55 

± 

0.60 

12 

LP  944-020 

18.69 

16.39 

13.98 

3 

10.725 

10.017 

9.548 

M9.0  V 

4 

0.08 

6.42 

± 

0.04 

7.04 

± 

1.32 

11 

GJ  1128 

12.74 

11.36 

9.62 

3 

7.953 

7.385 

7.037 

M4.0  V 

2 

0.15 

6.48 

± 

0.03 

6.33 

± 

1.00 

12 

DENIS  J 1048-3956 

17.37 

14.98 

12.47 

4 

9.538 

8.905 

8.447 

M8.0  V 

2 

0.08 

4.03 

± 

0.01 

4.48 

± 

0.73 

10 

SCR  1138-7721 

14.78 

13.20 

11.24 

4 

9.399 

8.890 

8.521 

M5.0  V 

3 

0.11 

8.31 

± 

0.07 

9.45 

± 

1.71 

12 

Proxima  Cen 

11.13 

9.45 

7.41 

3 

5.357 

4.835  4.384  M5.0  V 

Confirmed  Binaries 

2 

0.11 

1.30 

± 

0.01 

1.15 

± 

0.18 

12 

b 

LHS  1582AB 

14.69J 

13.33J 

11.60J 

4 

9.799J 

9.177J 

8.854J 

M4.5  VJ 

1 

19.67 

± 

0.47 

13.27 

± 

2.25 

12 

c 

GJ  748 AB 

11.10J 

9.95J 

8.47J 

3 

7.087J 

6.572J 

6.294J 

M3.5  VJ 

2 

10.00 

± 

0.12 

7.69 

± 

1.26 

12 

c 

LHS  3738AB 

15.78J 

14.29J 

12.46J 

3 

10.654J 

10.091J 

9.761J 

M4.5  VJ 

5 

19.15 

± 

0.38 

18.50 

± 

2.96 

12 

c 

Notes:  (a)  Distance  estimate  only  applicable  to  M  dwarfs;  (b)  Actual  error  is  ±  0.002  pc  in  trigonometric  distance;  (c)  “J”  signifies  joint  photometry  and  spectroscopy  for  unresolved  binaries. 
References.  For  spectral  types:  (1)  Reid  et  al.  (1995);  (2)  This  work;  (3)  Henry  et  al.  (2006);  (4)  Dieterich  et  al.  (2014);  (5)  Hawley  et  al.  (1996) 
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Figure  1.  Astrometric  residuals  plotted  in  right  ascension  (R.A.)  and  declination  (decl.),  in  units  of  milliarcseconds  (mas).  Filled  circles  represent  the  mean  of 
typically  three  to  five  frames  taken  in  a  single  night.  Open  circles  represent  nights  for  which  there  is  only  one  frame.  All  panels  are  on  a  ±50  mas  vertical  scale. 


The  goodness-of-fit  between  the  observed  data  yt  with  errors 
a,,  and  a  flat  line  with  y  0  was  determined  using  the  reduced 
chi  squared  statistic 


K 


±i 


(i) 


where  K  is  the  number  of  degrees  of  freedom,  given  by 
N  —  P  —  1  for  N  data  points  and  P  fitted  parameters.  Because 


we  did  not  attempt  to  fit  a  model  to  the  data,  but  only  to  analyze 
how  well  a  simulation  fits  the  data,  we  set  P  equal  to  zero.  As 
data  points  with  small  error  bars  are  more  heavily  weighted,  we 
discarded  epochs  that  have  unrealistic  errors  smaller  than  1  mas 
— representing  less  than  5%  of  the  epochs — because  such 
points  overconstrain  the  orbits  that  can  be  fit. 

Simulated  orbits  for  which  was  greater  than  4  in  at  least 
one  of  either  R.A.  or  decl.  we  consider  to  be  orbits  that  we 
would  have  detected.  This  threshold  of  4  is  based  empirically 
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Figure  2.  Periodograms  and  false  alarm  probabilities  (FAPs)  for  the  astrometric  residuals.  The  dashed  lines  correspond  to  FAP  values  of  0.1  (lowest  line),  0.01,  and 
0.001  (highest  line).  For  the  three  binaries,  the  periodograms  have  significant  peaks  (FAP  <  0.001)  at  the  same  period  in  both  R.A.  and  deck,  while  none  of  the  other 
targets  have  any  significant  peaks  even  at  the  0.01  level. 


on  the  values  plotted  in  Figure  3,  where  the  lower  panel  is  an 
inset  of  the  upper.  The  plotted  x  2  values  were  calculated  for 
(1)  the  stars  included  in  this  paper,  (2)  calibration  stars  with 
flat  residuals  we  use  to  monitor  potential  fluctuations  in 
equipment  and  our  data  reduction  pipeline,  and  (3)  additional 
known  binaries  with  perturbations.  We  compare  these  targets’ 
residuals  to  the  case  of  flat  residuals,  i.e.,  an  exact  astrometric 
solution  with  no  perturbation.  As  expected,  the  binaries  with 
perturbations  have  large  xr2d  values,  indicating  that  a  flat  line  is 

a  poor  fit  to  the  data.  The  two  solid  points  inside  the  \  2  =  4 
box  are  long  term  perturbations  with  large  gaps  in  the 


astrometric  observations.  As  more  data  are  collected,  those 
two  points  will  move  to  larger xL  values  outside  the  box,  but 
we  include  the  points  here  for  completeness.  For  those  stars 
with  no  perturbations,  all  have  xr2d  less  than  4  in  both  axes,  and 
are  centered  around  1.  This  indicates  that  we  are  accurately 
calculating  our  measurement  errors,  and  provides  an  empirical 
Wed  va^ue  °f  4,  above  which  we  are  sensitive  to  perturbations. 

4.  RESULTS 

Table  1  gives  the  parallax  and  proper  motion  results  for  the 
16  systems,  with  details  about  the  astrometric  observations 
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Figure  3.  Reduced  chi-squared  (\rJd)  values  comparing  the  astrometric 
residuals  to  the  case  of  perfectly  flat  residuals,  i.e.,  an  exact  astrometric  solution 
with  no  perturbation.  The  lower  panel  is  a  zoom  of  the  dotted  square  region  in 
the  upper.  Asterisks  represent  the  13  stars  without  close  stellar  secondaries 
analyzed  in  this  paper.  Diamonds  represent  additional  calibration  stars  with  flat 
residuals.  Solid  points  represent  binary  stars  with  perturbations,  including  the 
three  binaries  analyzed  in  this  paper,  which  are  circled.  All  of  the  targets  with 
confirmed  perturbations  have  values  greater  than  4  in  both  R.A.  and  deck, 
with  the  exception  of  two  stars  with  long-term  perturbations  that  have  gaps  in 
the  astrometric  observations.  All  the  targets  without  perturbations  have  x^d 
less  than  4  in  both  axes. 


(filters  used,  number  of  seasons  observed,  number  of  frames 
used  in  reductions,  time  coverage,  span  of  time,  and  the  number 
of  reference  stars)  and  results  (relative  parallaxes,  parallax 
corrections,  absolute  parallaxes,  proper  motions,  position 
angles  of  the  proper  motions,  and  the  derived  tangential 
velocities  based  on  relative  proper  motions  and  parallaxes).  All 
but  2  of  the  16  systems  have  parallax  errors  of  ~2  mas  or  less. 
BD-10  3166  and  GJ  876  have  larger  errors  due  to  combinations 
of  faint  reference  stars  and  short  exposures.  Corrections  to 
absolute  parallax  are  generally  less  than  ~2  mas,  so  systematics 
in  the  corrections  should  not  significantly  affect  the  results. 
Three  targets  (GJ  317,  GJ  667C,  and  BD-10  3166)  have 
corrections  of  ~  3^4  mas  due  to  reddening  of  the  reference 
stars,  which  skews  their  photometric  distance  estimates.  In 
these  cases,  we  adopt  a  generic  correction  of  1.50  ±  0.50  mas. 
The  per  observation  precision  for  each  target  is  listed  in 
Column  16,  representing  the  mean  of  the  observation  errors  in 
R.A.  and  decl.  The  percentage  of  companions  eliminated  listed 
in  Column  17  is  discussed  in  Section  5. 


Nine  of  the  16  targets  in  this  paper  have  parallaxes 
previously  published  by  RECONS,  and  are  noted  in  Column 
18.  The  results  presented  here  supersede  those  published 
previously  by  RECONS  because  additional  data  and  improved 
reduction  techniques  have  been  used,  as  discussed  in  detail  in 
Subasavage  et  al.  (2009) .  The  identical  parallax  of  LP  944-020 
is  also  presented  in  Dieterich  et  al.  (2014)  as  part  of  a  study  of 
the  stellar  hydrogen  burning  limit.  For  BD-10  3166,  we  did  not 
run  simulations  because  it  is  too  massive  and  far  away  for  us  to 
detect  any  type  of  substellar  companion.  However,  we  do 
provide  the  first  accurate  parallax,  and  conclude  that  BD-10 
3166  is  not  physically  related  to  the  star  with  a  similar  proper 
motion,  LP  731-076,  which  is  20"  away  (Bartlett  et  al.  2014,  in 
preparation). 

Photometric  and  spectroscopic  results  are  provided  in 
Table  2.  VRI  photometry  was  taken  using  the  CTIO  0.9  m 
(number  of  nights  of  photometry  in  Column  5),  with  errors  in 
VRI  <  0.03  mag  (Winters  et  al.  2011).  JHK  photometry  was 
retrieved  from  the  Two  Micron  All  Sky  Survey  (2MASS, 
Skrutskie  et  al.  2006)  catalog.  Spectral  types  are  given  in 
Column  9  with  references  in  Column  10.  Mass  estimates  were 
calculated  as  discussed  in  Section  3.  The  photometric  distances 
are  calculated  using  the  VRIJHK  distance  relations  (number  of 
relations  in  Column  14)  detailed  in  Henry  et  al.  (2004).  For 
systems  with  photometric  and  trigonometric  distances  that 
agree  within  the  errors,  we  conclude  that  they  lack  nearly  equal 
luminosity  companions.  The  trigonometric  distances  of  GJ 
748AB  and  LHS  1582AB  are  greater  than  their  photometric 
distances  due  to  companion  contributions  to  the  systems’  total 
flux.  The  two  distances  of  LHS  3738AB  agree  well,  indicating 
that  the  companion  is  significantly  fainter  than  the  primary.  GJ 
317  and  GJ  849  have  discrepant  (at  3.5  er  and  3.1  tr, 
respectively)  photometric  and  trigonometric  distances,  which 
does  not  necessarily  mean  that  these  stars  have  stellar 
companions,  as  main  sequence  stars  within  the  same  spectral 
type  can  vary  somewhat  in  luminosity. 

Bonfils  et  al.  (2013a)  note  a  radial  velocity  drift  in  their 
observations  of  GJ  849.  Montet  et  al.  (2014)  also  note  this 
drift,  and  constrain  the  minimum  companion  mass  to 
M  sin  i  <  2.5  MJup.  Our  astrometry  would  show  a  photocenter 
shift  for  unequal  mass  components,  as  discussed  in  Section  2. 
Only  components  of  roughly  equal  luminosity  and  mass  would 
provide  the  additional  flux  with  no  perturbation.  Such  a 
companion  would  have  been  observed  to  separations  as  close 
as  1"  in  our  images,  which  corresponds  to  ~9  AU.  At  a 
semimajor  axis  of  9  AU,  the  orbital  period  is  29.5  years  for  an 
equal  mass  companion.  This  results  in  velocities  for  each 
component  9.1kms_1  for  edge-orbits.  Thus,  for  most  orbital 
inclinations,  such  a  stellar  companion  is  ruled  out  by  the  radial 
velocity  data.  Therefore,  it  is  unlikely  that  a  stellar  companion 
similar  to  the  primary  is  contributing  to  the  overluminosity  we 
observe. 

Figure  4  shows  the  range  of  periods  and  masses  for  which 
90%  of  simulated  companions  would  have  been  detected,  based 
on  the  simulations  for  objects  with  masses  from  0.5  to  80  MJup. 
As  discussed  in  Section  3,  the  masses,  periods,  and  orbital 
parameters  of  the  simulated  companions  were  assigned 
randomly  from  a  uniform  distribution.  For  this  discussion  we 
set  the  dividing  line  between  planets  and  brown  dwarfs  at  13 
Mjup,  and  the  dividing  line  between  brown  dwarfs  and  stars  at 
80  Mjup.  The  bottom  panel  of  Figure  4  is  an  inset  of  the  top, 
showing  the  best  case  targets  in  more  detail.  The  noisy  nature 
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90%  Detection  Rates 


Figure  4.  Masses  and  periods  at  which  90%  of  companions  would  have  been 
detected,  based  on  simulations  of  10  million  companions  for  each  star. 
Companions  above  the  lines  would  have  been  detected.  The  lines  end  at  the 
length  of  observations  for  each  target.  The  six  M  dwarfs  with  known  planets 
are  labeled  in  the  top  panel.  The  bottom  panel  is  a  zoom  of  the  top,  showing  in 
greater  detail  the  six  additional  M  dwarfs  that  are  more  favorable  to  the 
astrometric  detection  of  planets.  GJ  1061  is  shown  as  a  dotted  line  to 
differentiate  it  from  GJ  1128.  Note  the  different  vertical  scales,  where  the  top 
panel  represents  primarily  brown  dwarfs  and  the  bottom  panel  Jovian  planets. 

of  the  lines  is  due  to  the  sizes  of  the  bins  used  in  the 
simulations.  The  bin  sizes  were  chosen  to  achieve  a  reasonably 
high  resolution,  while  still  having  enough  simulated  compa¬ 
nions  in  each  bin.  The  minimum  detectable  companion  mass  is 
smallest  for  companions  with  long  periods,  which  produce  the 
largest  amplitude  perturbations  in  the  astrometric  data.  The 
lower  panel  indicates  that  for  the  best  case  targets  (stars  at  close 
distances  and  of  low  mass).  We  are  most  sensitive  to  Jovian- 
type  planets  in  Jovian-type  orbits. 

For  a  companion  at  a  given  mass  and  orbital  period,  the 
amplitude  of  the  resulting  astrometric  perturbation  depends  on 
the  orbital  parameters  of  the  system,  and  the  mass  and  distance 
of  the  primary.  The  detection  limits  we  report  are  based  on 
simulations  of  companions  with  a  wide  range  of  masses, 
periods,  and  orbital  parameters.  Therefore  we  give  a  few 
representative  examples  of  how  the  results  in  Figure  4  translate 
into  astrometric  perturbations  in  mas.  In  the  case  of  a  face-on, 
circular  orbit,  a  20  MJup  companion  in  a  4  year  orbit  around  GJ 
581  would  cause  a  16  mas  perturbation,  while  a  15  Mjup 
companion  in  an  8  year  orbit  would  cause  a  20  mas 
perturbation,  and  a  10  Mjup  companion  in  a  12  year  orbit 
would  cause  a  17  mas  perturbation.  For  circular,  face-on  orbits 
around  Proxima  Centauri,  companions  of  1.5,  1,  and  0.5  MJup 


in  orbits  of  4,  8,  and  12  years  would  cause  perturbations  of  12, 
13,  and  8  mas,  respectively.  These  values  are  significantly 
greater  than  the  per  observation  precisions  listed  Table  1 — 
7.93  mas  for  GJ  581  and  4.83  mas  for  Proxima.  Thus,  the  90% 
detection  thresholds  given  in  Figure  4  are  reasonable 

For  the  four  planet  hosts  observed  longer  than  eight  years. 
Column  17  of  Table  1  gives  the  percentages  of  simulated 
brown  dwarf  companions,  ranging  from  81%-99%,  eliminated 
with  orbits  between  two  and  eight  years,  and  masses  between 
13  and  80  MJup.  Approximately  92%  of  all  simulated  brown 
dwarfs  have  been  eliminated  as  companions  to  those  stars 
known  to  host  exoplanets.  For  the  six  more  astrometrically 
favorable  targets,  we  calculate  the  percentages  of  simulated 
planetary  companions  eliminated  with  orbits  between  two  and 
eight  years  and  masses  between  1  and  13  MJup,  with  results 
ranging  from  69%-99%.  We  have  eliminated  ~86%  of  all 
simulated  planets  with  masses  of  1-13  MJup  around  these  six 
astrometrically  favorable  stars,  and  effectively  all  brown  dwarf 
companions  in  orbital  periods  of  2-8  years. 

Photocentric  orbital  solutions  for  the  three  binaries  are 
shown  in  Figure  5  with  the  corresponding  orbital  parameters 
given  in  Table  3.  From  our  astrometric  data  for  GJ  748  AB,  we 
find  an  orbital  period  of  2.504  ±  0.025  years,  which  is 
consistent  with  the  two  detailed  studies  of  the  system  by  Franz 
et  al.  (1999),  who  found  P  =  2.466  ±  0.008  years,  and 
Benedict  et  al.  (2001),  who  found  P  =  2.469  ±  0.001  years 
using  HST  Fine  Guidance  Sensor  data.  However,  we  determine 
an  eccentricity  of  0.06,  which  is  inconsistent  with  the  value  of 
0.45  found  in  both  of  the  HST  studies.  We  utilized  the  orbit¬ 
fitting  code  described  in  Hartkopf  et  al.  (1989)  and  set  starting 
eccentricities  of  0.05-0.95  in  increments  of  0.05;  regardless, 
our  data  converged  to  the  e  =  0.06  value  each  time. 

The  discrepancy  between  our  eccentricity  and  that  of  the 
HST  studies  is  likely  due  to  our  observations  of  GJ  748  AB 
having  been  taken  at  the  two  different  V  filters  discussed  in 
Section  2.1.  While  the  two  filters  are  photometrically  identical 
within  measurable  errors  (Jao  et  al.  2011),  they  are  not 
astrometrically  identical.  We  have  analyzed  the  astrometric 
residuals  for  over  500  targets  without  detectable  perturbations 
in  the  three  different  filters  ( V ,  R,  and  I )  over  the  length  of  our 
observing  program.  The  R  and  I  filters  are  stable,  but 
astrometric  offsets  in  the  V  filters  are  evident  over  the  time 
period  when  the  problematic  “new”  V  filter  was  used.  These 
offsets  have  been  mitigated  as  discussed  in  Section  2.1, 
allowing  data  from  both  V  filters  to  be  used  to  produce  reliable 
parallax  results.  In  the  case  of  GJ  748  AB,  we  are  able  to 
recover  the  correct  period,  but  the  offsets  in  the  residuals  are 
likely  contributing  to  the  errant  eccentricity.  We  presently  do 
not  possess  enough  observations  to  perform  a  reduction  of  GJ 
748  AB  without  the  problematic  V  filter  data.  We  include  the 
current  solution  because  it  is  our  only  system  for  which  an 
accurate  period  has  been  published,  against  which  to  compare 
our  results. 

In  contrast  to  our  photocenter  data,  the  FGS  observations 
resolve  the  system  into  two  components  at  15  epochs  over  1.8 
years.  They  are  exquisitely  sensitive  to  the  separation  and 
position  angle  of  the  secondary  from  the  primary,  and  are  to  be 
preferred  to  our  ground-based  results  for  the  eccentricity.  A 
clever  suggestion  by  Hugh  Harris  of  USNO  has  been  suggested 
to  solve  this  dilemma.  Because  in  an  unresolved  system  the 
center  of  mass  location  is  unknown,  the  zero  points  for  the 
residuals  in  R.A.  and  deck  are  unknown.  By  shifting  the  zero 
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Figure  5.  Top  panels  show  the  orbital  fits  for  the  astrometric  residuals  of  three  confirmed  binaries.  Filled  circles  represent  the  mean  of  typically  three  to  five  frames 
taken  in  a  single  night.  Open  circles  represent  nights  for  which  there  is  only  one  frame  and  are  not  used  to  derive  the  fit.  Bottom  panels  show  the  residuals  to  the  orbital 
fits.  All  panels  are  on  a  ±50  mas  scale. 


Table  3 

Orbital  Parameters  of  Known  Binaries 


Name 

P 

T0 

^  phot 

e 

i 

Long.  Peri,  (a;) 

Long.  Node  (0) 

(years) 

(mas) 

(deg) 

(deg) 

(deg) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

GJ  748 AB 

2.504  ±  0.025 

2005.86  ±  0.25 

28.1  ±  1.6 

0.06  ±  0.04 

137.8  ±  9.0 

218.2  ±  40.3 

173.6  ±  12.1 

LHS  1582AB 

5.309  ±  0.049 

2001.84  ±  0.14 

21.9  ±  1.3 

0.17  ±  0.03 

143.6  ±  7.5 

62.0  ±  14.3 

97.9  ±  12.4 

LHS  3738AB 

6.141  ±  0.059 

2005.73  ±0.16 

28.1  ±  1.2 

0.12  ±  0.02 

131.8  ±  4.1 

130.5  ±  11.2 

130.5  ±  5.1 

points  and  fitting  the  residuals,  a  different  eccentricity  may  be 
derived.  We  await  the  acquisition  of  resolved  data  for  several 
more  systems  before  exploring  this  technique  so  that  a  robust 
analysis  can  be  accomplished.  At  present  errors  on  the 
eccentricities  in  Table  3  should  be  treated  with  caution.  The 
orbital  solutions  for  LHS  1582AB  and  LHS  3738AB  are 
updated  and  improved  over  those  presented  in  Riedel  et  al. 
(2010),  which  were  the  first  orbits  presented  for  each  system. 
In  addition  to  demonstrating  the  astrometric  detection  and 
characterization  of  unresolved  companions,  these  results  can 
provide  additional  dynamically  determined  masses  for  M 
dwarfs,  once  the  systems  have  been  resolved. 

5.  DISCUSSION 

In  narrow-angle  field  astrometry  programs  such  as  the  one 
outlined  here,  the  lowest  mass  companions  detectable  by 
astrometry  are  dependent  on  a  number  of  factors.  These  include 
the  apparent  brightness  of  the  host  star  and  availability  of 
suitable  reference  stars,  which  affect  the  precision  of  the 


astrometric  measurements.  Additionally,  the  target  star  mass, 
companion  mass,  the  luminosity  ratio,  orbital  separation,  and 
system  distance  determine  the  size  of  the  astrometric  perturba¬ 
tion.  The  M  dwarfs  with  planets  that  have  been  discovered  so 
far,  including  the  six  in  this  paper,  are  not  the  most  favorable  to 
our  astrometric  observations.  They  are  among  the  brighter  red 
dwarfs,  with  the  latest  type  stars  at  M3. 5  V  for  the  radial 
velocity  detections  included  in  this  paper.  We  are  rarely 
sensitive  to  planets  around  these  stars.  However,  we  are  able  to 
rule  out  the  presence  of  a  large  fraction  of  potential  brown 
dwarf  companions  with  masses  of  13-80  MJup.  As  observations 
continue,  we  will  be  able  to  rule  out  companions  at  longer 
periods,  particularly  in  the  cases  of  GJ  317  and  GJ  1214,  which 
we  have  so  far  only  observed  for  3-4  years.  In  addition,  we 
have  demonstrated  that  we  are  sensitive  to  planets  with  masses 
of  1-13  Mjup  around  M  dwarfs  that  are  more  favorable  to 
astrometric  observations.  Of  particular  interest,  we  find  that 
Proxima  Centauri  has  no  companions  more  massive  than 
2  Mjup  with  periods  of  2-5  years  and  more  massive  than  1  Mjup 
for  5-12  years.  This  is  to  be  contrasted  to  the  results  of  other 
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Table  4 

Planets  Orbiting  M  Dwarfs  within  25  pc 


Name 

R.A.  (J2000.0) 

Decl.  (J2000.0) 

^mean 

N  n 

Ref. 

M  sin  (/) 

P 

Ref. 

Notes 

(mas) 

(M  Jup) 

(years) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

GJ  163  b 

04  09  15.6 

-53  22  25 

66.59  ±  1.79 

2 

1,2 

0.03 

0.024 

9 

GJ  163  c 

0.02 

0.070 

9 

GJ  163  d 

0.09 

1.654 

9 

GJ  176  b 

04  42  55.7 

+18  57  29 

110.00  ±  2.00 

2 

1,2 

0.03 

0.024 

10 

GJ  179  b 

04  52  05.7 

+06  28  35 

80.82  ±  3.78 

2 

1,2 

0.82 

6.264 

11 

GJ  317  b 

08  40  55.7 

-23  28  00 

65.34  ±  0.39 

3 

1,3,4 

1.81 

1.895 

4 

a 

LHS  2335  b 

10  58  35.0 

-31  08  38 

50.55  ±  1.55 

1 

5 

0.02 

0.007 

12 

GJ  433  b 

11  35  26.9 

-32  32  23 

112.09  ±  1.43 

2 

1,2 

0.02 

0.020 

13 

a 

GJ  1 148  b 

11  41  44.6 

+42  45  07 

88.81  ±  2.14 

2 

1,2 

0.30 

0.113 

14 

GJ  436  b 

11  42  11.0 

+26  42  23 

98.95  ±  2.07 

2 

1,2 

0.07 

0.007 

15 

GJ  581  b 

15  19  26.0 

-07  43  20 

159.28  ±  1.32 

3 

1,2,3 

0.05 

0.015 

16 

GJ  581  c 

0.02 

0.035 

16 

GJ  581  d 

0.02 

0.183 

16 

GJ  581  e 

0.01 

0.009 

16 

LP  804-027  b 

16  12  41.7 

-18  52  31 

69.46  ±  3.12 

1 

2 

2.10 

0.306 

17 

GJ  649  b 

16  58  08.8 

+25  44  39 

97.28  ±  1.32 

2 

1,2 

0.33 

1.638 

18 

a 

GJ  1214  b 

17  15  18.9 

+04  57  50 

69.04  ±  0.54 

3 

1,3,6 

0.02 

0.004 

19 

b 

GJ  667 C  b 

17  18  57.1 

-34  59  23 

138.24  ±  0.57 

3 

1,3,7 

0.02 

0.020 

13 

a 

GJ  667 C  c 

0.01 

0.077 

13 

GJ  674  b 

17  28  39.9 

-46  53  42 

220.11  ±  1.39 

2 

1,2 

0.03 

0.013 

20 

GJ  832  b 

21  33  33.9 

-49  00  32 

202.03  ±  1.00 

2 

1,2 

0.64 

9.353 

21 

GJ  849  b 

22  09  40.3 

-04  38  26 

113.46  ±  1.32 

3 

1,2,3 

0.90 

5.241 

22 

a 

GJ  876  b 

22  53  16.7 

-14  15  49 

214.45  ±  0.57 

4 

1,2, 3,8 

1.95 

0.167 

23 

GJ  876  c 

0.61 

0.082 

23 

GJ  876  d 

0.02 

0.005 

23 

GJ  876  e 

0.04 

0.342 

23 

Notes.  To  provide  a  uniform  format,  minimum  masses  and  orbital  periods  are  in  some  cases  rounded  to  fewer  significant  digits  than  in  the  original  publications,  (a) 
Additional  companion(s)  not  listed  on  both  exoplanet. eu  and  exoplanets.org;  (b)  transiting  planet;  all  others  are  radial  velocity  detections.  Recent  discoveries  GJ  191 
b/c  (Anglada-Escude  et  al.  2014)  and  GJ  687  b  (Burt  et  al.  2014)  were  not  listed  on  exoplanets.org  as  of  2014  July  1. 

References.  (1)  Yale  Parallax  Catalog,  van  Altena  et  al.  (1995);  (2)  Hipparcos,  van  Leeuwen  (2007);  (3)  This  work;  (4)  Anglada-Escude  et  al.  (2012);  (5)  Riedel 
et  al.  (2010);  (6)  Anglada-Escude  et  al.  (2013):  (7)  Fabricius  &  Makarov  (2000):  (8)  Benedict  et  al.  (2002);  (9)  Bonfils  et  al.  (2013);  (10)  Forveille  et  al.  (2009); 
(11)  Jioward  et  al.  (2010);  (12)  Bonfils  et  al.  (2011);  (13)  Delfosse  et  al.  (2013);  (14)  Haghighipour  et  al.  (2010);  (15)  Maness  et  al.  (2007);  (16)  Mayor  et  al. 
(2009);  (17)  Apps  et  al.  (2010);  (18)  Johnson  et  al.  (2010);  (19)  Harpspe  et  al.  (2013);  (20)  Bonfils  et  al.  (2007);  (21)  Bailey  et  al.  (2009);  (22)  Montetetal.  (2014); 
(23)  Rivera  et  al.  (2010). 


companion  searches  of  Proxima  Centauri.  Using  /  1ST  Fine 
Guidance  Sensor  data,  Benedict  et  al.  (1999)  report  a 
companion  detection  sensitivity  of  1  Mjup  at  a  60  day 
period.  Based  on  7  years  of  RV  observations,  Endl  & 
Kiirster  (2008)  preclude  the  presence  of  companions 
with  M  sin  (z)  ^  1  MNep  (0.05  Mjup)  at  periods  ^2.7  years. 
Together  these  studies  eliminate  all  Jupiter  mass  planets  around 
Proxima  Centauri  for  orbital  periods  out  to  12  years. 

In  the  broader  context,  these  results  are  consistent  with 
recently  published  searches  for  Jovian  companions  to  M  dwarfs 
at  shorter  orbital  periods.  Transit  searches  are  unlikely  to  detect 
companions  at  Jovian  orbits,  due  to  the  narrow  range  of 
detectable  inclinations.  Berta  et  al.  (2013)  found  no  Jupiter¬ 
sized  planets  in  their  transit  search,  which  is  most  sensitive  to 
companions  at  orbital  periods  less  than  10  days,  and  conclude 
that  such  planets  rarely  orbit  M  dwarfs.  Based  on  the  first  two 
years  of  their  astrometric  search,  Sahlmann  et  al.  (2014)  find  no 
planetary  mass  companions  to  the  20  M  and  L  dwarfs  they 
observed.  They  determine  the  occurrence  rate  of  planets  more 
massive  than  ~5  Mjup  to  have  an  upper  limit  of  9%. 

The  longer  time  coverage  of  the  astrometric  results  presented 
in  this  work  overlap  most  closely  with  radial  velocity  results. 


Bonfils  et  al.  (2013a)  report  detection  limits  based  on  radial 
velocity  measurements  for  102  nearby  M  dwarfs,  of  which  six 
are  featured  in  this  paper.  Among  the  102  M  dwarfs  searched, 
they  confirm  only  two  planets  with  orbital  periods  longer  than 
100  days.  At  an  orbital  period  of  1000  days  (2.7  years),  they 
report  a  detection  sensitivity  of  M  sin  (i)  <  2  MJup  around  90% 
of  stars  observed,  and  M  sin  (i)  <  30  Mjup  at  10,000  days 
(27  years).  Montet  et  al.  (2014)  report  that  6.5%  ±  3.0%  of  M 
dwarfs  host  a  1  to  13  MJup  planet  at  a  separation  less  than 
20  AU,  based  on  their  radial  velocity  survey  and  high 
resolution  imaging. 

The  results  presented  here  are  among  the  first  astrometric 
searches  for  Jovian  companions  at  Jovian  orbits,  and  fill  in 
relatively  unexplored  mass  and  period  parameter  space — 
models  of  how  M  dwarf  planetary  systems  form  and  evolve 
must  now  explain  the  lack  of  massive  companions  with  long 
period  orbits  like  those  in  our  Solar  System.  Looking  forward, 
the  Gaia  mission  may  detect  up  to  2600  planets  within  100  pc 
(Sozzetti  et  al.  2014).  However,  its  ability  to  determine 
accurate  masses  and  orbits  will  be  limited  to  orbital  periods  less 
than  6  years,  a  fraction  of  the  time  coverage  of  our  ground- 
based  astrometric  observations. 
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Finding  nearby  M  dwarfs  with  planets  remains  an  important 
challenge,  as  the  closest  planets  are  the  brightest  and  most 
easily  studied,  and  M  dwarfs  dominate  the  stellar  population. 
As  part  of  its  mission  to  characterize  the  solar  neighborhood, 
RECONS  is  developing  a  database  of  all  objects  with  accurate 
trigonometric  parallaxes  placing  them  within  25  pc  {nUl„  ^  40 
mas  with  an  error  ^10  mas).  For  an  extrasolar  planet  to  be 
included  in  the  RECONS  Database  it  must  orbit  a  star  that 
meets  the  above  criteria,  or  be  a  free-floating  object  with  a 
comparable  parallax,  and  be  listed  in  both  the  Extrasolar 
Planets  Encyclopaedia  (exoplanet.eu)  and  the  Exoplanet  Orbit 
Database  exoplanets.org,  Wright  et  al.  2011).  The  RECONS 
Database  currently  contains  1074  systems  having  M  dwarf 
primaries  (9.0^  My  <  21.0)  within  25  pc,  only  17  of  which 
have  detected  exoplanets  as  of  2014  July  1.  listed  in  Table  4. 
The  error-weighted  mean  parallax  for  each  system  is  given  in 
Column  4,  including  the  parallaxes  in  this  work  and  published 
values.  Minimum  masses  and  orbital  periods  for  planets  with 
references  are  listed  in  Columns  7-9.  Note  that  only  3  of  the  26 
reported  planets  have  masses  greater  than  Jupiter’s.  Based  on 
the  36  M  dwarf  primaries  within  5  pc  (Henry  2013),  we 
anticipate  that  there  are  4500  M  dwarf  primaries  within  25 
parsecs,  yet  only  17  so  far  have  been  found  to  host  planets. 
Clearly,  many  planets  lurk  undetected  in  the  solar  neighbor¬ 
hood.  Discovering  these  planets  will  require  a  wide  variety  of 
survey  techniques,  and  as  sensitivities  are  improved,  astrometry 
will  continue  to  play  an  important  role. 
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conversations  regarding  brown  dwarf  luminosities,  and  his 
assistance  in  surveying  the  existing  literature.  We  also  wish  to 
thank  Hugh  Harris  for  his  helpful  suggestion  regarding  the 
eccentricities  of  unresolved  binaries.  This  effort  has  been 
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efforts  of  the  National  Optical  Astronomy  Observatories  and 
the  members  of  the  SMARTS  Consortium.  This  research  has 
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Data  Explorer  at  exoplanets.org,  as  well  as  the  Extrasolar 
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